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Abstract 


The production of CO, gas at the DMFC anode leads to dramatic increases in pumping power requirements and reduced power output because 
of mass transfer limitations as bubble trains form in the channels of larger stacks. Experimental observations taken in a 5cm? DMFC test cell 
operated at 60°C, 1 atm, and with a methanol/water fuel flow rates of 5-10 cm? min“! indicate that the rate of bubble formation can be reduced 
by increasing the fuel flow because more liquid is available for the CO, to dissolve in. Further observations indicate that KOH and LiOH added to 
the fuel eliminates CO, gas formation in situ at low concentrations because of the greatly increased solubility that results. 

A mathematical model for the volumetric rate of CO2 gas production that includes effects of temperature and solubility is developed and extended 
to include the effects of hydroxide ions in solution. The model is used to predict the onset location of gas formation in the flow field as well as the 
void fraction at any point in the flow field. Predictions from the model agree very well with our experiments. Model predictions explain differences 
in the initial location of bubble formation for fuel solutions pre-saturated with CO, as opposed to CO>-free solutions. Experiments with KOH and 
LiOH added to fuel solutions confirm the validity of the model extension that includes solubility that is enhanced by chemical reaction. 

Experiments with LiOH, KOH, and ammonium hydroxide show that the long-term durability of standard Pt-Ru/Nafion®/Pt membrane electrode 
assemblies is compromised because of the presence of lithium, potassium, and ammonium cations that interact with the Nafion® membrane 
and result in increasing the ohmic limitations of the polymer electrolyte membrane. Experiments with Ca(OH),, while reducing gas formation, 


precipitate the product CaCO; out of solution too rapidly for downstream filtering, blocking channels in the flow field. 


© 2007 Elsevier B.V. All rights reserved. 
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1. Introduction 


Since its conception, the fuel cell has been attractive for the 
conversion of chemical energy into electrical energy due to high 
efficiencies and low emissions. While H2 powered cells have 
numerous potential applications, the need to power portable 
devices is driving development of fuel cells with liquid fuels. 

Liquid methanol has a higher energy density than hydrogen 
and other gas fuels and is more easily processed and stored. 
The direct methanol fuel cell (DMFC) takes advantage of these 
properties and is a potential power source for portable applica- 
tions; however, its current shortcomings in kinetics and methanol 
crossover problems yield a power density nearly an order of 
magnitude lower than that of a hydrogen fueled PEMFC [1]. 
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Because of those limitations, every opportunity to optimize cell 
operation must be considered. 

The DMFC consists of a polymer electrolyte membrane with 
anode and cathode catalyst layers on either side. Gas diffusion 
layers made of either carbon paper or carbon cloth are used to 
cover the catalyst layers to enhance diffusion and create elec- 
trical contact between the catalyst layers and current collectors 
on each end. The current collectors have some machined flow 
channel pattern in the surface that allows supply of fuel and 
oxidants. In a DMFC methanol is oxidized in the presence of 
water to CO2 at the anode while oxygen is reduced to water at 
the cathode. Operating temperatures are generally held below 
100°C to prevent boiling of the fuel. The following reactions 
describe the process: 


CH30H + H20 —> CO)+6H™ +6e7; 
(1) 


Anode reaction : 
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Cathode reaction: (3/2)O2+ 6H* + 6e” > 3H20; (2) 


Overall reaction : CH30H + (3/2)02 > CO2 +2H20. (3) 
To ensure maximum current density and optimal reaction condi- 
tions, efficient removal of both CO2 from the anode surface and 
water from the cathode surface is essential. While water man- 
agement at the cathode has been studied [2] due to its occurrence 
in Hy PEM cells, removal of CO3 at the anode has been largely 
overlooked and remains an important issue for research in the 
development of DMFCs. 

In the typical DMFC, methanol in a low concentration (1 M) 
solution is fed to the anode flow field and contacts the anode 
catalyst via the diffusion layer while the COz being produced 
traverses the diffusion layer to the anode flow field in reverse. 
The quantity of CO2 being produced is directly related to 
current produced by the cell. Combining Faraday’s constant 
with the anode stoichiometry from (1), CO2 in the amount of 
1.036 x 1074 mol min™! Amp~! is produced. The fuel for the 
reaction is liquid water and methanol, so the CO, produced at 
the anode is initially in the liquid phase. Because of the large 
quantities of CO that evolve, the liquid fuel quickly becomes 
saturated and gas phase CO2 forms. This gas formation, under 
poorly optimized operating conditions, may consume a signif- 
icant fraction of the flow field volume. CO, gas limits mass 
transfer of methanol and water to the anode catalyst surface and 
the presence of gas is indicative that the fuel solution is saturated 
with CO and that diffusion of CO2 from the catalyst surface may 
be limited. Initiation of two-phase flow leads to a capillary bub- 
ble train in the flow field creating pressure drops that decrease the 
efficiency of a stack by both increasing pumping power require- 
ments and methanol crossover. In portable systems the problem 
of fuel supply is generally overcome by increasing the flow rate 
of the anode fuel, which can have the undesirable side effects of 
further increasing methanol crossover and creating temperature 
gradients if the fuel solution is not properly pre-heated. Unfor- 
tunately, the formation of capillary bubble trains at high current 
production is inevitable even under these conditions. 

The topic of CO2 gas evolution in DMFCs has received 
some attention in the literature. Gas bubble growth and release 
from the diffusion layer was investigated by Mench et al. [3]. 
Argyropoulos et al. [4] witnessed the blocking of parallel flow 
channels at low flow rates and high current densities while study- 
ing the CO2 flow characteristics visually in situ. Yang et al. [5] 
conducted a comprehensive investigation into effects of flow 
field designs and cell orientation while qualitatively studying 
CO, bubble behavior using visualization in situ. While Yang et 
al. did not encounter the channel blocking slugs that Argyropou- 
los encountered, they did show that the presence of large gas 
slugs in the flow field correlated to decreased cell performance. 

In this paper, a model for the volume rate of CO2 gas produc- 
tion is presented which depends on the operating temperature 
of the cell, current produced by the cell, feed rate of fuel to 
the anode, and the solubility of CO% in the fuel as a function 
of temperature. Extensions of this model to include additional 
parameters that could contribute to gas reduction are briefly dis- 
cussed. Experimental results validating the original model are 


then presented as well as results from experiments exploring 
possible options for delaying the onset of gas formation. Reac- 
tants considered for gas onset delay enhancement include NaOH, 
LiOH, Ca(OH)2, and ammonium hydroxide. The long-term 
durability of standard Pt-Ru/Nafion®/Pt membrane electrode 
assemblies (MEAs) is investigated as it relates to these reactants. 

Under conditions when two-phase flow occurs, the ability to 
characterize it is important from a design standpoint so that prob- 
lems relating to heat transfer, mass transfer, and pressure drop 
all may be avoided. A key property of two-phase flow character- 
ization is the void fraction and, to this end, the model described 
in this paper is used to develop expressions determining exactly 
when gas will begin forming in the anode flow field. 


2. Theoretical expression for volumetric gas formation 
2.1. General expression 


The anticipated quantity of gas formation in can be related to 
the current discharge from the cell on a molar basis by 


dn izi 4 
dt 6F’ 
which is derived from Eq. (1) with the aid of Faraday’s constant, 
F, where n is the molar quantity of CO2 generated, and i, is the 
total output current divided by the length of the channel, L. 

If the liquid fuel is not saturated with CO2 prior to entering 
the cell, then there is some capacity for CO2 to remain dissolved 
in the liquid after its generation at the anode surface, thus not 
appearing in the gas phase. For the length of cell flow path where 
the fuel is not saturated with CO7, the appearance of bubbles will 
depend on the rate of CO2 formation by reaction and the local 
fluid mixing. If the CO2 production rate is sufficiently high and 
the bubbles are effectively nucleated, then CO2 gas phase will 
appear even if the fuel is not saturated with CO2. This would 
most likely occur at higher currents. 

In the region where all of the CO2 remains in the liquid phase, 
the change in concentration in the liquid can be calculated using 
a mass balance on a differential slice of the channel to get 


dCco, ix 


=—, 5 
fuel dx 6F (5) 


where fuel is the volumetric liquid flow rate, x the distance along 
the channel, and Cco, is the concentration of CO? in the liquid 
in mol1~!. Integrating Eq. (5) gives the CO2 concentration in 
the liquid as a function of distance along the flow path: 


iyX 


6Fqfuel 


Cco, (x) = (6) 
While high methanol concentrations can impact the solubility 
of CO, in ternary systems, a 1 M concentration typically used 
in a DMFC, the solubility can be approximated as water [6]. For 
the limited conditions where a typical DMFC is operated under 
atmospheric pressure in the temperature range of 298-373 K, 
Carroll et al. [7] gives detailed Henry’s law correlations for the 
solubility of COz in water. Considering Eq. (6) and the results 
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of Carroll et al., the distance into the channel where the liquid 
becomes saturated, xsat, will be 


6F teed CE, (T) 


lx 


(7) 


Xsat = 
where C&O, (D) is the saturation concentration of the liquid fuel 
solution. At distances greater than Xsat, gas will begin to form. 
Using the rate of CO2 production from Eq. (4) and the ideal 
gas law, the change in the volumetric gas flow rate, gco,, with 
respect to distance will be 


dqco, RT ix 

dx P 6F’ 
which can be integrated to give the gas flow rate at any location 
in the flow field channel: 


(8) 


_ $F atuei C&T) 


0 x< - 
lx 
4Co: = : i (9) 
i RT ix 6 Fafe Cto, (T) 
sk = Xsat) a a 
P OF lx 


Results from Eq. (9) are displayed in Fig. 1, where a temperature 
of 60°C, a pressure of 1 atm, and a current of 3 Amps over the 
entire flow path of arbitrary length, operated with various fuel 
flow rates. As the flow rate of the fuel is increased, the capacity 
to contain CO7 in the liquid phase gets larger, resulting in an 
extension in the delay of gas formation. 


2.2. Void fraction 


Because void fraction is a useful characterization of two- 
phase flow, taking Eq. (9) one step further results in 


aco l (10) 


ġ = gas fraction = | 
qfuel + 4C02 


which assumes that a negligible amount of liquid is consumed 
by the reaction. Fig. 2 displays the sudden rise in void fraction 


CO; gas flow rate, cc min“? 


- T T Ze T T 
0.0 0.2 0.4 0.6 0.8 1.0 
Fraction of flow distance 


Fig. 1. The gas flow rate predicted by Eq. (9) is displayed for a variety of fuel 
flow rates with respect to the fraction of the flow distance (T=60°C, P= 1 atm, 
i=3 Amps). The delay in the onset of formation of gas occurs because of the 
total solubility of CO2 in the fuel solution. 
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Fig. 2. The void fraction predicted by Eq. (10) is displayed for a range of fuel 
flow rates with respect to the fraction of the flow distance (T=60°C, P= 1 atm, 
i=3 Amps). Increasing the liquid flow rate is one way to reduce bubble forma- 
tion. 


after the initial onset of gas formation and is generated using 
the operating parameters of 60°C, | atm, 3 Amps, a flow field 
of arbitrary length, and various fuel feed rates. Under these con- 
ditions, small cells which should only require small amounts of 
fuel require high fuel flow rates simply to keep the void fraction 
in check. 


2.3. Delaying onset of CO2 gas formation 


The idea of the solution having a maximum carrying capac- 
ity for liquid-phase CO2 can be extended to include additional 
parameters having the potential for user manipulation. Despite 
possible detrimental effects on the methanol oxidation reaction, 
which occurs under acidic conditions, alkaline solutions are suit- 
able for illustration of how gas production might be reduced. In 
the case of the addition of a hydroxide salt, Eqs. (11) and (12) 
describe the additional chemical reactions [8]: 


CO2 + OHT = HCO3", (11) 
CO3” +H20 = HCO; +OH-, (12) 


where the equilibrium constants are Kj; =6.1 x 1071 mol! and 
Kı2=1.7 x 1074 mol1~!. At low concentrations an equilibrium 
is established between CO3*~ and HCO3~ obeying 


[HCO3"]__ 1.7 x 1074 
[C03] [OH] 


(13) 


Under conditions where OH is completely consumed, HCO37 
is the only species present, Eq. (13) is the only relevant reaction 
and OH” acts to solubilize CO» on a 1:1 mol basis. In this case, 
the concentration of OH” can be inserted into Eq. (7) to form 


6 Fqreea (Co, + Cop-) 


lx 


; (14) 


Xsat = 


where Coy is the molar concentration of hydroxide in the fuel 
solution. Here, the hydroxide concentration acts as an additional 
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sink into which CO? can exist in the liquid phase. Eq. (9) then 
becomes 


- 6 Fafe (CEO, (T) + Con-) 


0 x< A 
qco:= 4 - 
? RT ix 6Fqfuei (CEO, (T) + Con-) 
—— (x — Xsat) x > : 
P OF lx 


(15) 


Additional terms can be added to Eq. (7) or (14) in a similar 
fashion, as new methods for retention of CO7 in the liquid phase 
become available. 


3. Experimental 
3.1. Experimental setup 


The fuel cell hardware assembly used for the experiments 
is a 5 cm? model provided by Fuel Cell Technologies, Inc. The 
cell consisted of a pair of POCO graphite blocks with a pre- 
cision machined single serpentine flow-pattern and gold plated 
connectors fastened with aluminum end plates. Characteristics 
of the flow channel included a cross-sectional area of 1 mm?, 
an overall centerline length of 325 mm, an open ratio of 65.0%, 
and fourteen 180° turns with a diameter of 1 mm. 

Characteristics for membrane electrode assemblies used in 
all experiments included an active area of 5 cm?, an anode cat- 
alyst loading of 4.0 mg Pt-Rucm~”, a cathode catalyst loading 
of 2.0mgPtcm~?, and Nafion-117 polymer electrolyte mem- 
brane. The gas diffusion layers at the anode and cathode were 
carbon cloth and ETEK ELAT®, respectively. 

All experiments in this study used the setup shown in Fig. 3. 
Methanol solution isolated from atmosphere by a nitrogen pad 
was provided by a computer controlled GA series Micropump™ 
with steady flow rates available from 0.01 to 50cm? min~!. The 


DAQ & Control 
System 
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standard for all anode fuel solutions was HPLC grade methanol 
and deionized water (18 MQ). The methanol concentration for 
all experiments was 1 mol17!. 

Controlled extra dry oxygen without humidification was 
provided to the cathode by a feed system consisting of an 
MKS™ mass flow controller with flow rates available up to 
500 scc min~!. A purge gas consisting of ultra pure nitrogen was 
used to remove oxidant from the system prior to experiments. 

Gaseous CO, was separated and captured from unreacted 
methanol solution as it drained from the cell. Oxidant exiting 
the cathode was vented to atmosphere. 

Electrical conditions of the test cell were maintained by a 
Scribner 890CL fuel cell test load system controlled by a PC 
running Scribner’s proprietary FuelCell® (Version 3.4d). All 
tests reported here were performed using a cathode feed rate 
of 100 scc min~!, the cell in a vertical orientation, a current dis- 
charge control mode, and atmospheric pressure (101.325 kPa). 


3.2. Model verification procedure 


Quantification of the rate of CO2 gas generation for model 
verification was accomplished through the use of a gas capture 
device. CO2 gas separated from the unspent fuel exiting the 
anode was passed into a chamber isolated from atmosphere. 
Separating the chamber from atmosphere was a pad of water pre- 
saturated with CO> to prevent additional unintended absorption 
of CO, gas. The rate of gas generation was then determined 
by measuring the volume displacement of water with respect to 
time. Data were taken in sets at various current discharge rates 
while keeping the operating temperatures, pressures, and fuel 
flow rates constant. 

For verification of the model extension, where hydroxide salts 
were added to the standard anode fuel solution, the cell was 
operated with a standard anode fuel solution until it stabilized, 


Scribner 890CL 
Test Load 


§ Methanol 
È Solution 
SEH 2 
MicroPump 
i S] ] i 
DMFC CO,(g) 
H,0() + 0,(9) | COA) * | Gas-Liquid 
CH,OH(/) + H,O0(/) O i 
Water CH,OH(/) 
+ H,O(/) 
N, O, 


Fig. 3. Schematic of the DMFC test loop and CO7 gas capture system. 
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whereupon the source was switched to modified fuel and the 
above procedure followed. New membrane electrode assemblies 
were used for both LiOH and KOH experiments. 


3.3. Stability experiments with hydroxide salt and amine 
solutions 


Knowing that the addition of an alkaline compound might 
interfere with the methanol oxidation reaction and that it might 
adversely affect the Nafion® membrane, experiments were con- 
ducted to determine the stability of the membrane electrode 
assemblies under conditions where hydroxide salts were present 
after model verification experiments showed promising results. 
NaOH, LiOH, Ca(OH)2, and ammonium hydroxide were each 
tested in separate experiments using fresh membrane electrode 
assemblies. Each test was a total of 65min long and was 
conducted at a temperature of 60°C, an anode flow rate of 
5cm? min™!, and at a current discharge rate of 400 mA cm~? 
(2 Amps). Each test began with standard fuel solution, cur- 
rent was stepped to 400 mA cm~? from open circuit, 5 min was 
allowed to pass for equilibrium to occur, and fuel was switched 
to a solution that included an addition of hydroxide salt that 
resulted in a salt concentration of 0.005 M. After the fuel switch, 
the electrical performance of the cell was monitored for 60 min. 


4. Results and discussion 


Eq. (9) predicts a stepwise function where the rate of gas 
formation increases at a rate proportional to current, but is not 
present until some onset current where the formation of CO2 
exceeds the carrying capacity of the fuel solution. The carrying 
capacity of the fuel solution is the maximum rate at which the 
fuel solution can absorb CO2, which is a function of the fuel 
flow rate and the solubility of CO? in the solution, which itself 
is a strong function of both temperature and pressure. While 
applying the ideal gas law to Eq. (4) alone would predict an 
increase in gas formation for a temperature rise of 10°C due 
to Charles’ law, it would not consider the effect of flow rate. 
Despite a competing effect of temperature on solubility, the flow 
rate of the fuel solution, if appropriately raised, can compensate 
for the increase in temperature and actually decrease the rate of 
gas formation. Fig. 4 compares the results of model verification 
experiments with predictions of both Eq. (9) and the ideal gas 
law applied to Eq. (4) while simultaneously demonstrating the 
shift produced by the delayed onset of gas formation and the 
overwhelming effect of fuel flow rate. In practice, small bubbles 
are occasionally witnessed earlier in the flow field than the model 
predicts (Yang et al.). This is a result of localized CO? saturation 
within the gas diffusion layer. 

Yang et al. states that an increase in the sweeping rate of 
bubbles due to an increase in flow rate causes a reduction in the 
number of bubbles in the flow field. This can be misinterpreted 
to mean that that an increase in flow rate causes a reduction in the 
volume consumed by gas, which Fig. 4 rejects. An assumption 
of steady state requires that the quantity of CO7 exiting the cell 
must equal what is generated, which is dependant solely on the 
cell’s current discharge rate. Without consideration for the CO» 
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Fig. 4. The delay of gas formation due to CO? solubility is illustrated by plot- 


ting the rate of CO» generation as a function of total current (@: T=60°C, 
Q=5 cm? min7!, pre-saturated fuel solution; ©: T=60°C, Q=5 cm? min7!; 
v: T=70°C, Q=10cm? min“ !; (---) Eq. (5) prediction; (—) Eq. (9) model 


predictions). 


solubility of the fuel solution, this entire quantity of CO2 must 
immediately evolve as gas, which, using the ideal gas equation 
as the sole correction factor, would produce a greater volume of 
gas to form regardless of the flow rate and contrary to the data. 

Adjustments to temperature and flow rate may be prohibited 
by external requirements of the cell, or they may not be suffi- 
cient to completely eliminate gas under high current discharge 
conditions. Additionally, changes to flow rate and temperature 
may cause more problems due to pumping requirements and 
methanol crossover than they solve by the reduction of two- 
phase flow, thus is is of interest to enhance the solubility of CO2 
in the fuel system. 

Eq. (15) predicts that additional control over the effective 
solubility of CO3 in the fuel solution can be achieved by the addi- 
tion of components that react with CO2. Of particular interest 
are reactions involving hydroxides and amines, both well under- 
stood and used in industrial CO2 stripping operations. Potassium 
hydroxide, tested in a concentration of 0.0053 M (Fig. 5), causes 
a shift in the rate of gas formation that matches the prediction 
of Eq. (15). It is worth pointing out that the addition of KOH in 
a concentration of 0.0053 M produces a nearly identical shift to 
that of increasing the flow rate in Fig. 4. Similar reduction occurs 
with LiOH, which is again predicted by Eq. (15). Increasing the 
concentration to 0.01 M, the rate of gas production continues to 
decrease, and along with it, the onset of gas formation is delayed 
further (Fig. 6). 

Despite significant reduction in CO2 formation and model 
verification, when hydroxide was added the cell did not remain 
stable. Experiments were conducted to investigate the stability of 
membrane electrode assemblies due to the previous observations 
and to explore other reaction possibilities. After an initial drop 
in power density at 2 Amps, the cell using ammonium hydrox- 
ide experienced a linear decay of 0.45 mW cm~? min~! over the 
final 35 min of the experiment, resulting in an overall 34% drop 
in power density. Sodium hydroxide has a similar rate of decay 
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Fig. 5. The delay of gas formation when KOH is added to the fuel solution is 
illustrated by comparing the rate of CO2 gas generation as a function of cur- 
rent (@: T=60°C, Q=5 cm? min7!, pre-saturated fuel solution; ©: T=60°C, 
Q=5cm? min“!; Y: T=60°C, Q=5cm? min™!, 0.0053 M KOH; (- - -) Eq. (5) 
prediction; (—) Eq. (15) predictions). 


at 0.37 mW cm~? min! over the final 35 min, however it expe- 


riences a larger initial drop, resulting in a 33% overall drop in 
power density for the experiment. Lithium hydroxide fares much 
better with a decay rate nearly half that experienced by ammo- 
nium hydroxide at 0.23 mW cm~? min™! which, with inclusion 
of the initial drop, results in a 23% decrease in power density. 
Calcium hydroxide, exhibits somewhat different behavior: there 
is an initial rise in current density, followed by a steep decline 
with instabilities growing in magnitude. Fig. 7 compares the 
results of these studies. Interactions of the hydroxide salt cations 
with the Nafion® membrane, which has proton conducting sites 
that can potentially undergo permanent replaced by +1 cations, 
is the likely explanation for the reduction in power output. 
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Fig. 6. The delay of gas formation when LiOH is added to the fuel solution is 
illustrated by comparing the rate of CO2 gas generation as a function of cur- 
rent (@: T=60°C, Q=5 cm? min7!, pre-saturated fuel solution; OQ: T=70°C, 
Q=10cm? min~!; Y: T=70°C, Q= 10cm? min“ !, 0.01 M LiOH; (- - -) Eq. (5) 
prediction; (—) Eq. (15) predictions). 
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Fig. 7. The power density for cells operated with selected hydroxide salts is plot- 


ted against time for identical operating conditions of T= 60°C, Q =5 cm? min™!, 


i=2 Amps, and [OH~]=0.005 M. 


Compelling visual evidence in Fig. 7 indicates that the results 
of the calcium hydroxide experiment differed from sodium, 
ammonia, and lithium counterparts. The calcium cation was not 
expected to cause the same permanent replacement of proton 
sites in the Nafion® membrane due to the +2 state of the calcium 
ions. The result of the reaction is insoluble calcium carbonate, 
which immediately precipitates out of solution eventually block- 
ing the flow field and any other part of the fuel path downstream 
of the cell. The rapid rate of reaction and, consequently, the rapid 
rate of solid formation, prevent effective downstream filtering 
and are the cause of the instabilities visible in Fig. 7. 


5. Conclusions 


A model that explains the rate of CO2 gas formation at the 
anode of a DMFC including a function of the cell’s output cur- 
rent, operating temperature, operating pressure, fuel flow rate, 
and the solubility of CO% in the fuel solution, which is also a 
strong function of temperature is presented and verified. The 
onset delay of gas formation is a function of solubility, which 
may be exploited by DMFC designers to reduce or possibly 
eliminate the formation of gas by properly tuning the operating 
parameters of the cell. Additionally, the expression developed 
to describe the rate of gas formation can be extended to include 
additional parameters capable of manipulation for design pur- 
poses. 

For the first time, the solubility of CO2 in DMFC fuel is 
increased by addition of OH. The results, again in agreement 
with the model, show that CO2 gas can be reduced or poten- 
tially eliminated through the addition of components to the fuel 
system. 

Long-term durability of the membrane electrode assembly 
was used as screening criterion in an exploration of possible 
candidates for extending the delay of gaseous CO? onset. While 
hydroxide salts were promising candidates due to their perfor- 
mance in model validation experiments, declines in power output 
of over 20% in an hour are unacceptable for any benefit received. 
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Encapsulating the ions potentially harmful to the MEA in such 
a way such that reactions are still allowed to take place, possibly 
with liquid membrane emulsions, might be a good avenue of 
future study. 
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